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FLAME CHEMILUMINESCENCE RATE CONSTANTS FOR 
QUANTITATIVE MICROGRAVITY COMBUSTION DIAGNOSTICS 

Jorge Luque,* Gregory P. Smith, Jay B. Jeffries/ David R. Crosley 

Molecular Physics Laboratory, SRI International 
Menlo Park CA 94025 

Absolute excited state concentrations of OH(A), CH(A), and C 2 (d) were determined in three low 
pressure premixed methane-air flames. Two dimensional images of chemiluminescence from 
these states were recorded by a filtered CCD camera, processed by Abel inversion, and calibrated 
against Rayleigh scattering. Using a previously validated 1-D flame model with known 
chemistry and excited state quenching rate constants, rate constants are extracted for the reactions 
CH + 0 2 — » OH(A) + CO and C 2 H + O -> CH(A) + CO at flame temperatures. Variations of 
flame emission intensities with stoichiometry agree well with model predictions. 


INTRODUCTION 

Optical emission signatures from combustion 
experiments provide a convenient chemically based 
diagnostic for a variety of applications. For example, 
microgravity experiments can be monitored without the 
instrumental complications of laser techniques. 
However, successful implementation of an emission 
diagnostic requires a model-based interpretation, and 
quantitative results depend on improved knowledge of 
the identity and rate constants for the reactions 
producing the excited electronic states observed. 

We report here experiments in low pressure premixed 
flames to determine these rate constants for CH(A), 
OH(A), and C 2 (d) production during methane 
combustion. A camera with appropriate optical filters 
records spatially resolved two-dimensional profiles of 
the emissions. An Abel inversion converts these 
images into flame-center height profiles. Rayleigh 
scattering from nitrogen is recorded with the same 
system to calibrate the measurements. Using a 1-D 
flame model which has been validated by flame radical 
measurements, previous LIF determinations of CH(A) 
and OH(A) quenching rate constants, and these results, 
we then derive rate constants for the reactions CH + 0 2 
— > OH(A) + CO and C 2 H + O — * CH(A) + CO. Three 
25-30 torr methane-air premixed flames of different 
stoichiometry are examined in this first work. 
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Modem laser diagnostics methods are typically used to 
provide species-specific, spatially-resolved 
measurements in combustion experiments, with the 
goal of providing some detailed data concerning the 
local progress of the flame chemistry, or the state of the 
fluid dynamics that has delivered a detected degree of 
combustion. 1,2 PLIF images and local Raman 
measurements are examples. But because sophisticated 
optical and laser equipment, which are often in need of 
adjustment or repair, require high electrical power, and 
are quite massive to transport and accommodate, are 
required to conduct such measurements, simpler 
experimental diagnostics would find good use in 
micro gravity and space-based flame experimentation. 
Optical flame emissions, like the familiar blue color of 
a gas stove, provide such an opportunity. This light is 
the result of specific molecular electronic transitions 
from excited states produced in non-equilibrium 
concentrations by specific chemical reactions. Thus the 
resulting chemiluminescence is a marker for this 
particular chemistry, and only simple optical detection 
equipment is required. These advantages were 
recognized early in the development of combustion 
science, and efforts were undertaken to characterize the 
spectroscopy and the reactions producing identified 
emissions. Gaydon provides an excellent summary of 
the emission spectra, flames, and possible reactions. 3 
For hydrocarbon-air flames, most light comes from 
OH(A-X) at 280-310 nm, CH(A-X) at 430 nm, CH(B- 
X) at 390 nm, Cjfd-a) Swan bands found mostly 
between 470-550 nm, the C0 2 flame band continuum 
between 320-500 nm, and polynuclear aromatic 
molecules and soot luminescence in rich systems 
emitting in the visible. 
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Little quantitative work on flame chemiluminescence 
has been published since Gaydon’s review, but 
microgravity combustion experimental programs, 4 and 
possible uses in combustion control applications 5 have 
increased interest recently. As one example of a 
chemical correlation, Najm et al 6 suggested the OH(A- 
X) ultraviolet emission measures the final stages of the 
CH, reduction chain, for example. Walsh et al 7,8 
quantitatively measured the optical emission intensities 
from CH(A) at 431 nm and OH(A) at 308 nm in a 
methane jet diffusion flame, and modeled the results. 
Similar experiments are being undertaken for these 
flames in microgravity environments. 4 In addition, they 
measured temperature and major species using Raman 
and Rayleigh scattering, and ground state CH, OH, and 
NO with laser induced fluorescence. 

The reliability of the model comparison, and its 
validation for use on the flight experiments, depends 
largely on the availability of reliable kinetics for the 
production of the excited state species that emit. The 
loss rates through collisional quenching are largely well 
known from recent flow cell, flame, laser pyrolysis, and 
shock tube experiments and theoretical formulations. 9 
The current model of Walsh’s data 7-8 errs in predicting 
the results by factors of 2-6, so quantitative success and 
interpretation for these diagnostics depends on 
obtaining improved measurements of the CH(A) and 
OH(A) production reactions. 

Studies of hydrocarbon flame chemiluminescence 
kinetics since Gaydon’s book are limited, and vary in 
their recommendations. Flame studies are presented by 
Joklik et al 10 and Grebe and Homann measure emission 
in a flow system. 11 Shock tube measurements are 
presented by Hwang et al, 12 and Mertens. 13 The only 
direct rate measurements of reactions producing CH(A) 
are by Devriendt et al 14 at room temperature. Among 
the older flame studies, suggested chemiluminescence 
rate constant values are given by Porter et al for OH 
and CH, 15 and by Bulewicz et al for CH. 16 

Low pressure premixed laminar flames provide an ideal 
setting to determine these production rate constants. A 
series of flames investigated in our laboratory have 
been well characterized by temperature and absolute 
OH, CH, NO, and HCO concentration distributions 
versus height above the burner surface, 17 and our 
understanding of these flames and their chemistry 
validated by agreement with one-dimensional flame 
calculation results using the GRI-Mech 3.0 
mechanism. 18 Quenching rates for excited states were 
measured by LIF, 9 and optical calibration techniques 
for absolute concentration determinations are 


operational. 19 In the study reported here, we measure 
absolute emission intensities and excited state 
concentrations for OH(A), CH(A), and C 2 (d) in three 
low pressure methane-air flames of varying 
stoichiometry. Using modeling results, we will then 
provide suggested rate parameters for production 
reactions for the OH and CH chemiluminescence. 


EXPERIMENTAL 

The three premixed methane-air flames studied here are 
those labeled Standard (<(>=1.07, 25 torr), Lean (<{>=0.81 , 
25 torr), and Rich (<|>=1.28, 30 torr) in previous LIF 
studies of methane flame chemistry. 17 Details, 
temperatures, and CH profiles versus height above the 
burner may be found there. 17 Operation at low pressure 
allows good spatial resolution of the species, by which 
the progress of the flame chemistry can be followed. 
The flames are supported on a 6 cm diameter sintered 
plug McKenna burner, housed inside a vacuum 
chamber designed for optical probing. A concentric 
flow of argon is used to isolate and stabilize the flame. 
The burner translates vertically with respect to the 
optical axis of the laser and detection system. A servo 
control valve maintains the chamber and flame 
pressure, while mass flow controllers regulate the 
methane, oxygen, and nitrogen. A laser beam traverses 
the chamber through Brewster angle windows, with 
power monitored at the exit (Laser Precision Rj-7200), 
for the Rayleigh scattering calibration of the optical 
detection system. A large silica window, through which 
light may be imaged onto either the spectrometer or 
CCD camera detection systems, lies perpendicular to 
the laser path. 

In order to compare experimental results for the excited 
state species concentrations and derive production rate 
constants, we performed one-dimensional flame model 
calculations using the Sandia code Premix 20 and the 
GRI-Mech 3.0 chemical mechanism for methane 
oxidation. 18 Measured flow rates and temperature 
profiles are also inputs. This follows the same 
procedure used for LIF measurements of ground state 
species in these flames. Reactions are added for the 
new excited state species OH(A) and CH(A) - radiative 
decay, collisional quenching, 9 and potential production 
steps. In addition, a few reactions involving C and C 2 
species were added; since these species may be 
involved in the production kinetics, a more complete 
description is desirable. These addenda to GRI-Mech 
3.0 are shown in Table 1, where the * signifies the 
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excited state. The excited state production rates are 
provisional and will be discussed next. 

The likely source of OH chemiluminescence cited by 
Gaydon 3 * and others is the reaction CH + 0 2 . (Other 
energetically feasible reactions are CH 2 or l CH 2 + 0 2 . 
Several H0 2 reactions and O + HCO can be ruled out 
by the observation that OH(A) is found at higher flame 
temperatures.) Porter et al 15 measured species and 
emissions in several 18 torr methane and acetylene 
flames, and using appropriate transition strengths and 
estimated quenching rates deduced the rate constant 
adopted for our mechanism. In hydrogen flames and in 
burnt gas regions, the recombination or inverse 
predissociation reaction O + H + M is responsible for 
OH(A) formation. Both Kaskan and Davis et al. 22 have 
derived rate constants from hydrogen flame data, 
assuming the reaction to be H + 20H — > OH(A) + H 2 0 
(a recombination assisted excitation). A reanalysis of 
some of their data in terms of an O + H + M reaction, 
using current quenching rates, 9 provides the tabulated 
rate constant, with an uncertainty of at least 60%. 
Experiments with our low pressure hydrogen flames 25 
could refine this value in the future. 

The reaction responsible for CH chemiluminescence, 
the A-X at 430 nm and B-X at 390 nm, is less certain. 
Three reactions have been proposed, and literature 
analysis of experiments often conflict or omit the other 
possible steps. For the reaction C 2 H + 0 2 , we made an 
Arrhenius fit to the room temperature rate constant of 
Devriendt et al. 14 and the shock tube emission results at 
2000-2600K from Hwang et al. 12 This is 4 times higher 
than Merten’s recent 3000K shock tube results, 13 which 
would suggest a higher activation energy and as much 
as a 16-fold reduction at 1500K flame temperatures. 
For the reaction C 2 H + O, at room temperature, 
Devriendt et al. 1 4 measure 1.1 x 10 13 cn^mole'V 1 but 
Grebe and Homann 11 propose a much lower value 
(dependent on their uncertain production kinetics for 
C 2 H) from their flow reactor data, 7.0 x 10 M , Joklik et 
al 10 estimate a value of 7.0 x 10 12 from their low 
pressure acetylene flame study. An average value of 
6.2 x 10 12 is used in the preliminary mechanism, with 
no temperature dependence assumed and large error 
limits likely at flame temperatures. Finally, the step C 2 
+ OH is favored as the source of the excited CH 
observed in flames by Bulewicz et al. 16 and Porter et 
al. 15 The mechanism rate constant is from a reanalysis 
of Porter’s results using a shorter revised C 2 (d) 
radiative lifetime 26 and our quenching rates for CH(A), 9 
and is about twice the estimated value given by 
Bulewicz. 16 


Modeling work on the source reaction rates for C 2 (d-a) 
emission must wait for measurements on the quenching 
loss rates, using LIF. Appearance of these emissions in 
rich flames or regions suggests to many that possible O 
atom reactions are not likely to be responsible. The 
isotope work of Ferguson 27 has shown that reactions 


Table 1 . Chemiluminescence Kinetics Additions 
k = AT"exp(-E/RT) A(cm 3 , mole,s) E(kcal/mole) 
(el2 = x 10 12 ) 

* T 1 C 1 1 


Reaction 

A 

n 

E 

Source 

C2+H2=C2H+H 

4.00e5 

2.4 

1000 

est. 

(C2H+H2) 

CH+CH=C2+H2 

5.00el2 

0 

0 

est. 

C+C+M=C2+M 

3.00el4 

0 

-1000 

est.(2X+M) 

C+CH=C2+H 

S.00el3 

0 

0 

est.(0+CH) 

0+C2=C+CO 

5.00el3 

0 

0 

est.(0+CH) 

C2+02=C0+C0 

9.00el2 

0 

980 

ref. 21 

CH+02=CO+OH* 

6.00el0 

0 

0 

ref. 15 

C2+OH=CO+CH* 

1.1 lei 3 

0 

0 

ref. 15 
reanalyzed 

C2H+0=CO+CH* 

6.20el2 

0 

0 

ref. 10,1 1,14 

C2H+02=C02+CH* 

4.10el3 

0 

4500 

ref. 12, 14 

0+H+M=OH*+M 

3.63el3 

0 

0 

ref. 22 
reanalyzed 

C+H+M=CH*+M 

3.63c 13 

0 

0 

est. 

(O+H+M) 

CH*-»CH 

1.86e6 

0 

0 

ref. 23 

CH*+N2=CH+N2 

303 

3.4 

-381 

ref. 9 

CH*+02=CH+02 

2.48e6 

2.14 

-1720 

ref. 9 

CH*+H20=CH+H20 

5.30el3 

0 

0 

ref. 9 

CH*+H2=CH+H2 

1.47el4 

0 

1361 

ref. 9 

CH*+C02=CH+C02 

0.241 

4.3 

-1694 

ref. 9 

CH*+CO=CH+CO 

2.44el2 

0.5 

0 

ref. 9 

CH*+CH4=CH+CH4 

1.73el3 

0 

167 

ref. 9 

OH*— »OH 

1.45c6 

0 

0 

ref. 24 

OH*+N2=OH+N2 

1.08ell 

0.5 

-1238 

ref. 9 

0H*+02=0H+02 

2.10el2 

0.5 

-482 

ref. 9 

0H*+H20=0H+H20 

5.92el2 

0.5 

-861 

ref. 9 

OH*+H2=OH+H2 

2.95el2 

0.5 

■444 

ref. 9 

0H*+CO2=OH+CO2 

2.75el2 

0.5 

-968 

ref. 9 

OH*+CO=OH4CO 

3.23el2 

0.5 

-787 

ref. 9 

OH*+CH4=OH+CH4 

3.36el2 

0.5 

-635 

ref. 9 


of two separate single-carbon-atom species are 
responsible. Candidates include C + CH, CH + CH, 
and C + CH 2 . It will be very hard to distinguish these 
experimentally, and may not matter operationally. 
Emission and LIF measurements in hydrogen arcjets 
seeded with methane may provide additional useful 
data. 28 

Low resolution sample spectra from the center region of 
maximum luminescence were recorded in the three 
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flames, using a 1.05 m focal length (f/4.5) ultraviolet 
camera lens to image the light onto the fiber optic 
entrance channel of a spectrum analyzer (SC 
Technology Process Sentry 400, resolution 3.2 nm). 
The results are shown in Figure 1 , after correcting for 
detector wavelength response. The main features are 
clearly identifiable as OH(A-X), CH(B-X), CH(A-X), 
and C 2 (d-a) bands. The CH and especially C 2 emissions 
are very dependent on stoichiometry, and the C 2 Swan 
bands cannot be identified in the lean flame spectrum 
(see the lOx insert). Some underlying emission, 
particularly near 350 nm and to the red of 430 nm, 
especially for the lean flame, may be due to the C0 2 * 
continuum. But it is too faint for worthwhile analysis. 
In the following analysis, we assume light from 280- 
340 nm is from OH(A-X), that between 380 nm and 
450 nm is excited CH, and beyond 450 nm for the 
standard and rich flames only is the C 2 Swan bands. 
Note the Av = +2 Swan bands will contribute to the 
intensity assigned to CH(A-X) near 430 nm, but little 
error is introduced - the contribution is probably less 
than 10% of the Av = +1 peak at 470 nm. 29 

A set of 2-D flame emission images were recorded 
through the same camera lens from 1 m(camera f/5.6) 
to obtain line-of-sight integrated luminescence signals 
for the three emitters and three flames. Signal was 
detected with a gated intensified CCD camera operating 
in simple shutter mode(Princeton Instruments, ICCD 
576G.RBT, 14 bits dynamic range, 384x576 pixels). 
Imaging an illuminated piece of graph paper established 
the pixel size to be 0.0135 cm at the flame, and 
identified the burner center and surface in the images. 
Five filters were used in front of the camera: a Schott 
GG455 glass long wave pass (LWP) filter for the Swan 
bands; a Schott GG375 LWP filter which will transmit 
CH as well as C 2 emissions; alternatively, a narrow 
band 430 ± 10 nm interference filter for CH(A-X); a 
narrow band 300 ± 10 nm interference filter passing 
some of the OH(A-X) light; and a UG-11 Schott glass 
filter which transmits light readily from 275-370nm to 
view all the OH(A-X). The difference image between 
the GG375 and GG455 filters covers the 370-460 nm 
window to isolate the CH emissions. Filter 
transmission curves were generated from spectra taken 
using the Process Sentry 400 spectrum analyzer and the 
outer cone excited OH and C0 2 uv-blue light from a 
Bunsen flame. The OH emission feature in this flame 
lets one determine the fraction of OH(A-X) radiation 
transmitted by the narrow band filter. The intensifier 
gain was calibrated from repeat measurements with 
different settings, and the manufacturers' wavelength 
dependence was assumed. This contributes about a 
factor of 3 difference when relating emitting OH to CH, 


and adds uncertainty to the OH(A) determination. (The 
optical system was calibrated at 434 nm .) Finally, in 
addition to these gain and transmission terms, the 
recorded emission intensities must be multiplied by the 
ratio of observation time to radiative lifetime in order to 
obtain excited state concentrations. If we have the 
same amounts of steady state OH(A) and CH(A) in a 
flame, we will see more CH emission because its 
radiative rate is faster. 

The images measure integrated emissivities through the 
depth dimension of the cylindrical flames, but 
comparison to 1-D chemical models requires localized 
excited state concentrations near the uniform center of 
the flame, as a function of height above the burner. An 
Abel inversion of the data must be performed, and we 
used the code of Walsh et al. 7 based on work by 
Dasch. 30 Binning of data pixels and operation at lower 
spatial resolution is necessary to avoid excessive noise 
in the inverted result. The current analysis added 16 
horizontal and 2 vertical pixels, and used 3 point 
smoothing, yielding a resolution of 0.05 cm height and 
0.43 cm radial. After inversion, the results provide 
emission intensities for each species and flame from 
specific local volume elements. 

To relate the point emission intensities of the Abel 
inverted images to excited state concentrations, the 
system optical detection efficiency must be calibrated. 
The same Rayleigh scattering method previously 
employed to determine absolute ground state 
concentrations from flame LIF measurements 19 was 
used. The CCD camera system with GG-375 filter is 
operated with a 50 ns gate triggered by the firing of an 
excimer pumped dye laser tuned to 434 nm. The laser 
beam is directed through the chamber above the burner, 
filled at 50-100 torr nitrogen, and the power (20 
pj/pulse at 10 Hz) is measured(Laser Precision Rj-4200 
meter). About 600 shots are averaged, and then the 
process is repeated with the nitrogen and laser beam 
replaced by the standard flame and the CH and C 2 
luminescence. The Rayleigh signal is integrated over a 
2.16 cm (160 pixel) length and the full 6 pixel height- 
width of the laser beam. 

The information is available to compute the number of 
Rayleigh excitations in the observation volume, N*(R), 
relate it to the observed signal intensity, and use this 
ratio to determine the CH(A) concentration from the 
intensity of CH(A-X) signal. 

N*(R) = 4no p d I(R) , (1) 
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where a = 2.4 x 10' 27 cm 2 is the N 2 Rayleigh cross- 
section at 434 nm, 1 d = 2.16 cm, p at 100 torr is 2.17 x 
10 18 cm 3 , and I(R) = 4.3 x 10 13 photons/pulse. (The 
camera is perpendicular to the laser, which is polarized 
in the plane of the optical system.) Then the excited 
state concentration 

[CH*] =s N(CH*)/V(Abel) 

= N*(R) x S(Abel)/S(R) x t(rad)/t(obs) (2) 

where the ratio of radiative lifetime to observation gate 
t(rad)/t(obs) = 537 ns / 50 ns, S are the respective signal 
levels, V(Abel) is the volume of the Abel inverted cell, 
N(CH*) is the number of excited state CH molecules in 
the cell, and N*(R) is from (1). 

REm T S 

Figure 2 displays the image of the OH(310 nm) light 
emission from the standard 25 torr 0=1.07 methane-air 
flame. The resulting Abel inversion is shown below the 
image. Results for the CH(370-460 nm) and C 2 (X>460 
nm) emissions are very similar. A flattening of the 
profiles can be seen, as might be expected from the 
inversion and the general saucer shape of species 
distributions in these flames. 31 Vertical cuts through 
these figures to compare height distributions show that 
the inversion has minor effects on most excited state 
peak positions, widths, and amounts. If one 

assumes, for example, that the excited CH is uniformly 
distributed over the nominal 6 cm burner diameter 
instead of taking the radial data and performing the 
inversion, only a 5% error is introduced. 

Traces of an uneven pattern are apparent in some of the 
Abel results. The halves are slightly asymmetric, the 
central region 2-4 cells often have 10-15% lower 
concentrations than the peak found toward mid-radius 
of the flame, and sometimes the central region is 
abruptly shifted toward the burner surface. For OH*, 
optical depth calculations suggest that absorption of 
emission from the flame center by OH(X) may 
contribute to this pattern. Future analysis refinements 
may reduce this -10% uncertainty from the inversion. 

The results for the maximum excited state 
concentrations and ratios between flames are given in 
Table 2, along with comparisons to the starting model 
results. The reaction C 2 H + 0 2 dominates model CH* 
production and this model greatly overpredicts the 
experimental results. A second model result without 
this step is given in the penultimate column, and gives 
better predictions. The temperatures at the 
chemiluminescence peaks range from 1500-1950K; at 


1700K and 25 torr, a concentration of 10 7 cm' 3 
corresponds to 70 parts per trillion. We estimate 
uncertainties of 20% in relative values, and 35% in 
absolute concentrations. The CH(A) interference filter 
results are consistent with the values obtained from the 
difference image of the glass filters shown in Table 2. 
The CH(A)/CH* ratios show 80% of the emission 
coming from the A state (0.8 for the standard and rich 
flames, 0.6 for the lean one), which roughly agrees with 
the Fig. 1 spectra. 


Table 2. Maximum Excited State Concentrations (cm* 3 ) 
i ( e8 = x 10 8 ) 


Species 

Flame(s) 

Exper. 

Model 

w/o 

C2H+02 

Expert 

Model 

OH* 

Rich 

2.33e8 

9.03e7 


2.58 

OH* 

Standard 

1.56e8 

5.38e7 


2.90 

OH* 

Lean 

1.20e8 

3.79e7 


3.17 

CH* 

Rich 

5.57e7 

2.1 le9 

4.4 le7 

.026,1.26 

CH* 

Standard 

9.62e6 

3.09e8 

9.93e6 

.031,0.97 

CH* 

Lean 

2.36e6 

4.66e7 

1.62e6 

.051,1.46 

C2* 

Rich 

1.06e7 




C2* 

Standard 

1.26e6 




OH* 

Rich/ 

Standard 

1.49 

1.68 



OH* 

Lean/ 

Standard 

0.77 

0.70 



CH* 

Rich/ 

Standard 

5.8 

6.8 

4.4 


CH* 

' 

Lean/ 

Standard 

0.24 
(-17 A) 

0.15 

i 

0.16 



The flame ratios of excited state concentrations are 
fairly well reproduced by the model results. This 
agreement supports the chosen chemistry and affirms 
the ability to predict the stoichiometry dependence of 
the emissions. The lean CH ratio is high for the 370- 
460nm glass filter result, but matches predictions for 
the narrow band CH(A-X) value, suggesting that some 
non-CH background may contribute (see also the Fig. 1 
spectrum). The rich CH* ratio lies between the values 
from models dominated by C 2 H + 0 2 and C 2 H + O, 
which suggests both reactions may contribute. With 
respect to absolute concentrations, the model 
underpredicts OH(A) by a factor of 3, and overpredicts 
CH* by up to a factor of 40. If the reaction C 2 H + 0 2 
— » CH* + C0 2 is removed, however, good agreement is 
obtained. 

The Abel inversions also provide height profiles of the 
excited state concentrations at flame center for 
comparison with model results. The results for OH* 
and CH* in the 3 flames are shown in the first two rows 
of Figure 3. Concentration scales are 10 s and 10 7 cm' 3 
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respectively, and decrease with the fuel fraction. The 
model curves are plotted for 3 times the OH* 
predictions, 1/30 the CH* predictions, and for the CH* 
formed from the O + C 2 H reaction only. 

The zero height location of the profiles derived from 
the images is uncertain. Repeat runs can differ by 0.05 
cm, and the average peak position is unexpectedly 0.1 
cm above the model prediction. Previous agreement of 
the model with CH(X) LIF peak positions, 17 the OH(A) 
precursor, implies the model should predict the 
chemiluminescence peak positions well. Some data is 
available that places the peak of the CH* emission with 
respect to the CH LIF maximum in the standard 
flame, 32 and it gives a height of 0.60 ± .02 cm for the 
maximum CH* emissivity. The nominal peak of the 
current CH* emission data (before inversion) is 0.72 
cm. Therefore all profiles shown in Fig. 3 have been 
adjusted backward by 0.12 cm. 

The last row of Fig. 3 shows the results for C 2 * in the 
rich and standard flames. The comparison model results 
- the product of C and CH concentrations - reflect one 
likely production reaction. 

mission 

The general order of chemiluminescence peaks in these 
flames, for both observation and model, is OH* before 
C 2 * and then CH*. The model widths appear slightly 
wider than the Abel-inverted measurements. Two 
model predictions are included for the CH* results, 
corresponding to C 2 H + 0 2 and C 2 H + O production 
reactions. The C 2 H + O runs provide a better fit for 
peak position, and amount. In general, the adjusted 
standard flame profiles show good agreement with the 
model, while the rich flame model is up to 0.08 cm later 
(higher) in the flame than the data. The lean flame 
OH(A) peaks later than the model, but a second 
experimental run showed much less disagreement. 

The interference filter results, labeled CH(A), are 
plotted in the Lean CH* panel along with the glass filter 
370-460 nm profile. The tail in the broadband profile at 
heights into the burnt gases, as well as the larger than 
expected difference in peak magnitudes (mentioned 
earlier), indicates a significant minority of the 
broadband emission (-20%) is not associated with 
CH*. This interpretation also improves model 
agreement, in this panel of Fig. 3 and for the Lean CH* 
rows of Table 2.. 

The starting model predicts about 30 times too much 
excited CH, and the production rate constants would 


need to be reduced by this fraction to obtain agreement 
(see Fig. 3). Analysis of the model shows that the C 2 H 
+ 0 2 reaction is responsible for nearly all CH* 
production. Even reducing this rate by 16, as Mertens’ 
work indicates, 13 is insufficient. A second model 
without this reaction predicts the measured amounts. 
(The multiplier for this model value is 1.1 ±0.2 .) 
Here, the C 2 + OH reaction only contributes about 10%; 
for this to be a significant cause of the 
chemiluminescence, the rate constant would need to be 
10 times higher and nearly gas kinetic - an unlikely 
proposition. Our conclusion is to recommend the rate 
constant of 6.2 x 10 12 cn^mol'V for O + C 2 H as the 
CH* source, although any combination of O and 0 2 
reactions with the same net rate will account for our 
results. Determining this division will require further 
experiments in flames with significantly different 0/0 2 
ratios. The rich/standard flame CH* ratio results in 
Table 2 would suggest nearly equivalent participation. 
Note that the O atom rate constant value is nearly the 
same as Joklik’s acetylene flame number (7 x 10 12 ), and 
indicates a factor of 2 decline from the measured room 
temperature value. 14 

The Walsh diffusion flame results 7,8 for CH* are 
overpredicted by a factor of 3 by a model that uses 
GRI-Mech 2.11 chemistry and the chemiluminescence 
rate constants from Devriendt et al. 14 Although there 
are enough other modeling uncertainties that one cannot 
derive precise values from this data, it would suggest an 
O + C 2 H rate constant of 3.6 x 10 12 , only 40% below 
our recommendation. 

The results for OH*, a multiplier of 2.9 ± 0.3, suggest a 
faster rate constant for CH *f 0 2 by a factor of 3, to 1.8 
x 10 u cm 3 /mole/s. Uncertainties are at least 35%, and 
the OH* result depends on the relative gain of the CCD 
camera in the ultraviolet (currently taken from the 
manufacturer’s specifications). Modeling of the Walsh 
et al diffusion flame OH* observations 7,8 using GRI- 
Mech 2.11 and the same OH* mechanism implies a 
multiplier of 1.6; the predicted position of the 
chemiluminescence is poor, however. The starting 
mechanism rate constant is from Porter’s 15 analysis of 
18 torr methane and acetylene flames; it depends on the 
ratio 1 cq[M][OH*]/[CH][ 0 2 ] , where k Q is the rate 
constant for quenching by M, [0 2 ] is determined by 
mass spectrometric sampling, and [OH*] and [CH] 
come from emission intensity and absorption 
measurements and respectively depend inversely on the 
transition probabilities. So k kg / [0 2 ] F(CH)F(OH) . 
The Porter analysis assumed an average quenching 
cross section of 7 A 2 , while our estimates with more 
recent quenching kinetics 9 range from 3.5 A 2 for the 
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acetylene flames to 9.0 A 2 for methane flames similar 
to our own. 9 Grebe and Homann 11 determined a 298K 
rate constant of 4.8 x 10 10 in an H/0/C 2 H 2 low pressure 
flow system, only 20% below Porter et al. 15 

The OH* concentrations and production rate constant 
depend on the camera intensifier gain relative to CH*. 
A direct calibration will refine the current preliminary 
values. The height profiles of Fig. 3 are also uncertain. 
Point measurements of emissivities and LIF as the 
burner height is scanned will provide greater accuracy. 

Future experiments in other low pressure flames, 
including C 2 H X fuels, will refine these rate constant 
values. Flames with N 2 0 oxidizer and hydrogen flames 
seeded with hydrocarbons will provide the varied 0/0 2 
ratios needed to decipher the production reaction(s) for 
CH*. Experiments with C-2 and C-3 fuels will extend 
the applicability of quantitative chemiluminescence 
diagnostics to higher hydrocarbon flames. 

LIF quenching measurements of C 2 (d-a) in low pressure 
flames will provide the decay information needed in the 
mechanism to deduce excited state C 2 production rate 
parameters from the emission data. In addition, C 2 (a) 
populations can also be determined to provide a test of 
the kinetics mechanism involving precursors to CH*. 

CONC LUSIONS 

Absolute excited state concentration profiles for 
OH(A), CH(A and B), and C 2 (d) were determined from 
optical emission intensities for low pressure methane- 
air flames of different stoichiometries. Production rate 
constants were then obtained from computer modeling 
of the flames using known excited state quenching 
rates. Stoichiometry variation is well predicted. The 
reactions and preliminary rate constants are: 

CH + 0 2 — » OH* + CO k = 1.8 x 10 n cm 3 /mole/s; 

O + C 2 H — » CH* + CO k = 6.2 x 10 12 cm 3 /mole/s. 

The results do not, however, rule out the C 2 H + 0 2 
reaction as an alternative source of CH*. These rate 
parameters can be used in modeling chemiluminescence 
measurements in methane flames to derive quantitative 
chemical information from experiments. 
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Figure 1 . Flame emission spectra from the center of the three low pressure premixed methane-air flames 



Figure 2. Bottom: 2-D Image of OH(A-X) emission from the d> = 1.27 methane-air 25 torr flames; Top: Abel 
inversion results giving the OH(A) concentration distribution through flame center as a function of radius and 
height. Bands show intensity contours to the light maxima at the flame center. Black background is < 10%. 


Figure 3. Experimental excited state profiles for OH(A), CH(A.B), and C 2 (d) (top to bottom) at the center of three 
low pressure methane-air flames of <t> = 1 .27. 1 .07. and 0.81 (left to right). Four modified modeling results as 
described in the text are shown for comparison. The experimental profiles have been shifted 0.12 cm closer to the 
nominal burner surface (see text). 
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